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The tremendous pandemic potential of coronaviruses was 
demonstrated twice in the past few decades by two global outbreaks 
of deadly pneumonia. Entry of coronaviruses into cells is mediated 
by the transmembrane spike glycoprotein S, which forms a trimer 
carrying receptor-binding and membrane fusion functions’. 
S also contains the principal antigenic determinants and is the 
target of neutralizing antibodies. Here we present the structure 
of a mouse coronavirus S trimer ectodomain determined at 4.0A 
resolution by single particle cryo-electron microscopy. It reveals 
the metastable pre-fusion architecture of S and highlights key 
interactions stabilizing it. The structure shares a common core with 
paramyxovirus F proteins~”, implicating mechanistic similarities 
and an evolutionary connection between these viral fusion proteins. 
The accessibility of the highly conserved fusion peptide at the 
periphery of the trimer indicates potential vaccinology strategies 
to elicit broadly neutralizing antibodies against coronaviruses. 
Finally, comparison with crystal structures of human coronavirus 
S domains allows rationalization of the molecular basis for species 
specificity based on the use of spatially contiguous but distinct 
domains. 

Coronaviruses are enveloped viruses responsible for 30% of mild 
respiratory infections and atypical pneumonia in humans worldwide’. 
The emergence of the severe acute respiratory syndrome coronavirus 
(SARS-CoV) in 2002 and of the Middle East respiratory syndrome 
coronavirus (MERS-CoV) in 2012 demonstrated that these zoonotic 
viruses can transmit to humans from various animal species, and sug- 
gested that additional emergence events are likely to occur. The fatality 
rate of SARS-CoV and MERS-CoV infections are about 10-37%'* and 
there are no approved antiviral treatments or vaccines. 

Coronaviruses use S homotrimers to promote cell attachment and 
fusion of the viral and host membranes. S determines host range, cell 
tropism and is the main target of neutralizing antibodies during infec- 
tion’. S is a class I viral fusion protein synthesized as a single chain 
precursor of about 1,300 amino acids that trimerizes upon folding. 
It is composed of an amino-terminal S,; subunit, containing the 
receptor-binding domain, and a carboxy-terminal S, subunit, driv- 
ing membrane fusion. Cleavage by furin-like host proteases at the 
junction between Sj and S) (S2 cleavage site) occurs during biogenesis 
for some coronaviruses such as mouse hepatitis virus (MHV, the 
prototypical and best-studied coronavirus)!°. The S; and S, subunits 
remain non-covalently associated in the metastable pre-fusion 
S trimer. After virion uptake by target cells, a second cleavage is mediated 
by endo-lysosomal proteases (S,’ cleavage site), allowing fusion activation 
of coronavirus S proteins®. 

Crystal structures of coronavirus S post-fusion cores demonstrated 
that the fusogenic conformational changes lead to the formation of 
a so-called trimer of hairpins that is the hallmark of class I fusion 
proteins’ '°. These structures contain two heptad-repeat (HR) regions 


present in S, assembled as an extended triple helical coiled-coil motif 
(HR1) surrounded by three shorter helices (HR2). Crystal structures of 
several coronavirus S receptor-binding domains in complex with their 
cognate receptors have also been reported'’’. Finally, cryo-electron 
microscopy (cryoEM) of SARS-CoV virions provided a snapshot of 
the S glycoprotein at 16 A resolution!». The lack of high-resolution data 
for any coronavirus S trimer has prevented a detailed analysis of the 
infection mechanisms. 

We produced an MHV S ectodomain trimer with enhanced stability 
by mutating the S, cleavage site and fusing a GCN4 trimerization motif 
at the C-terminal end of the construct. The resulting MHV S ecto- 
domain forms a trimer binding with high-affinity to the soluble mouse 
CEACAM 1a receptor (Extended Data Fig. 1a, b). We used state-of-the 
art cryoEM'® to determine the structure of the MHV S ectodomain 
trimer at 4.0 A resolution (Fig. la—c and Extended Data Figs 2 and 3). 
We fitted the crystal structures of two S; domains'’!*’” and built 
de novo the rest of the polypeptide chain using Coot'® and Rosetta!?”° 
(Fig. 1d-f, Extended Data Figs 2-4 and Supplementary Tables 1 and 2). 
The final model includes residues 15 to 1118, with an internal break 
corresponding to a loop immediately upstream from the S,’ cleavage 
site (residues 827-863). The region connecting the S, and S, subunits 
(residues 718-754) features weak density that correlates with its accessi- 
bility for proteolytic cleavage in vivo. Residues 453-535 were modelled 
by density-guided homology modelling using Rosetta owing to the 
poor quality of the density in this region (Extended Data Fig. 3k). 

The MHV S ectodomain is a 140 A long trimer with a triangular 
cross-section varying in diameter from 70 A, at the membrane proximal 
base, to 140 A at the membrane distal end (Fig. 1d, e). The structure 
comprises two functional subunits (Fig. 2a—d): a distal moiety con- 
stituted by the S; subunits; and a central stem connecting to the viral 
membrane formed by the S, subunits. 

The S; subunit has a “V’ shape contributing to the overall triangular 
appearance of the S trimer (Extended Data Fig. 5a). The S; N-terminal 
moiety comprises domain A, which is folded as a galectin-like B-sandwich 
decorated with extended loops on the viral membrane distal side, and 
a three-stranded antiparallel 3-sheet plus an a-helix on the viral mem- 
brane proximal side. The S; C-terminal half folds as three spatially 
distinct 8-rich domains, termed B, C and D (Fig. 2a-d). 

The S, subunit connects to the viral membrane and is characterized 
by the presence of long a-helices (Figs 2b-d and 3a). A central helix 
(C39) stretches 75 A along the three-fold molecular axis towards the 
viral membrane (Fig. 3a). It is located immediately downstream of the 
HRI1 motif, which folds as four consecutive a-helices (a2¢6-C9; Fig. 3a 
and Extended Data Fig. 6a, b), in sharp contrast to the 120-A long 
HRI helix observed in the post-fusion S structures’~’ (Extended Data 
Fig. 6c—-e). The 55-A-long upstream helix (029), so named because it 
is located immediately upstream of the S,’ cleavage site, runs parallel 
to and is zipped against the central helix via hydrophobic contacts 
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largely following a heptad-repeat pattern. A core antiparallel 3-sheet 
(346-B49-B50) is present at the viral membrane proximal end and is 
assembled from an N-terminal 3-strand (846), preceding the upstream 
helix, and a C-terminal 8-hairpin (849-650), located downstream of 
the central helix. 

MHV S, features a topology similar to the paramyxovirus F proteins 
(such as respiratory syncytial virus (RSV) F: root mean squared devi- 
ation (r.m.s.d.) 4A over 125 residues), with a comparable 3D organ- 
ization of the core 3-sheet, the upstream helix and the central helix 
(Fig. 3a, b). Importantly, these motifs were shown to remain invariant 
in the pre- and post-fusion F structures’. The conservation of these 
motifs among coronavirus S and paramyxovirus F proteins suggests that 
these fusion proteins have evolved from a distant common ancestor. 
Although the density is too weak to trace the polypeptide chain down- 
stream from 359, secondary structure predictions suggest that the domain 
directly preceding HR2 could adopt a similar fold in coronavirus S 
and paramyxovirus F proteins. 

In the S trimer, the three central helices are packed via their cen- 
tral portions whereas the two ends splay away from the three-fold 
axis (Extended Data Fig. 7a—c). Additional contacts between the 


Figure 1 | 3D reconstruction of the MHV S trimer determined 
by single-particle cryoEM. a-c, 3D map filtered at 4.0 A 
resolution coloured by protomer. Two different views of the 

S trimer (from the side (a) and from the top, looking towards the 
viral membrane (b)), and a side view of one S protomer (c) are 
shown. d-f, Ribbon diagrams showing the MHV S atomic model 
oriented as in a-c. 


upstream and central helices participate to inter-protomer interac- 
tions. Furthermore, the S; subunits interlock to form a crown around 
the S, trimer stabilizing it in the pre-fusion conformation (Fig. 3c, d 
and Supplementary Table 3). This is illustrated by the large surface 
area buried at the interface between each S, subunit and the S subu- 
nits of the three protomers (1,970 A”). Many of these contacts involve 
the HR1 helices and the fusion peptide region. These polypeptide seg- 
ments undergo major refolding during the fusogenic conformational 
changes (Extended Data Fig. 6a—e), which supports the notion that 
the S; subunits maintain the S, fusion machinery in its metastable 
state. Substitutions of the conserved alanine 994 by valine in helix 
Q2g or of the conserved leucine 1062 by phenylalanine in the cen- 
tral helix were shown to attenuate fusogenicity*’*. Our structure 
suggests that the former substitution would strengthen hydrophobic 
packing against the core B-sheet (Extended Data Fig. 7b), and that 
the later substitution could reinforce molecular stapling of the cen- 
tral helices (Extended Data Fig. 7a, c). The expected modification 
of the energy landscape between pre-fusion and post-fusion con- 
formations would explain the reduction in fusion activity of these 
mutants”), 


a S, subunit S, subunit Figure 2 | Architecture of the MHV S protomer. 

1036- a, Schematic diagram of the S glycoprotein 
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helix; CT, cytoplasmic tail; FP, fusion peptide; 
HR1/HR2, heptad-repeats; TM, transmembrane 
domain; UH, upstream helix. b-d, Ribbon 
diagrams depicting three views of the S protomer 
coloured as in a. Asterisk denotes the MHV S 
receptor-binding region. Disulfide bonds are 
shown as green sticks except for residues 453-535, 
for which they are not shown. 
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The predicted fusion peptide includes the C-terminal half of helix 
21 and extends up to the N-terminal half of «22 (refs 6 and 23) (Fig. 2c). 
Q1 is an amphipathic helix located at the periphery of the S trimer, 
burying hydrophobic side chains towards the S, centre and exposing 
charged residues to solvent (Fig. 2c and Extended Data Fig. 7b, c). In 
the case of porcine epidemic diarrhoea coronavirus, trypsin processing 
at the S,’ site can only occur after host cell attachment‘. This indicates 
that receptor binding could allosterically increase the accessibility of 
the S,’ site, which is located within helix a2). The acidic pH of the endo- 
lysosomes could also contribute to exposing the S,’ cleavage site for 
coronaviruses requiring cleavage in this compartment. The fact that 
helix a2; appears dynamic and is found immediately downstream from 
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Figure 3 | Pre-fusion structure of the coronavirus 
fusion machinery. a, b, Topology and ribbon 
diagrams showing the structural similarity 
between coronavirus MHV S, (starting at residue 
755) (a) and paramyxovirus RSV F (PDB 5C6B) 
(b). For clarity, only part of RSV F is shown, with 
conserved secondary structural elements coloured 
identically as for MHV Sp. ‘#’ denotes motifs 
participating to the post-fusion HR1 coiled-coil. 
c, d, Two different views of the MHV S trimer 
(from the side (c) and top, looking towards the 
host cell membrane (d)) highlighting how S, 
(ribbon diagram and semi-transparent surface) 
wraps around the S, fusion machinery (ribbon 
diagram) to stabilize it. 


a disordered loop suggests that it could undergo considerable ‘breathing’ 
motions. Regardless of the mechanism promoting cleavage, the MHV S 
structure reported here explains the requirement for processing at the 
S,’ site, as it frees the fusion peptide from the S, N-terminal region, 
which is a prerequisite for its insertion ~200 A away in the target mem- 
brane. The peripheral position of the fusion peptide is similar to what 
has been observed in the parainfluenza virus 5 F* and HIV gp41 (ref. 25) 
prefusion structures (Extended Data Fig. 8a—c). The notable accessibility 
of the fusion peptide and its sequence conservation among corona- 
viruses®”? suggest that it would be an ideal target for epitope-focused 
vaccinology initiatives aimed at raising broadly neutralizing antibodies 
against S glycoproteins (Fig. 4a—c and Extended Data Fig. 9). Major 
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Figure 4 | Potential strategy for neutralizing coronavirus infections. 
a, Surface representation of the MHV S trimer coloured according to 
sequence conservation using the alignment presented in Extended 
Data Fig. 9. The fusion peptide sequence is highly conserved among 
coronavirus S proteins. b, Surface representation of the MHV S trimer 


highlighting the peripheral position of the fusion peptide (blue and 
cyan). c, Ribbon diagrams of the MHV S trimer showing the overlapping 
positions of the fusion peptide (residues 870-887, blue and cyan) and of a 
major antigenic determinant identified for MHV and SARS-CoV (residues 
875-905, magenta spheres). 

00 MONTH 2016 


| VOL 000 | NATURE | 3 


© 2016 Macmillan Publishers Limited. All rights reserved 


RESEARCH §aSaia 


antigenic determinants (inducing neutralizing antibodies) of MHV 
and SARS-CoV S proteins overlap with the fusion peptide region and 
support the suitability of this approach”®”’. Antibodies binding to this 
site will not only hinder insertion of the fusion peptide into the target 
membrane, but will also putatively prevent fusogenic conformational 
changes. This epitope-focused strategy has proven successful to obtain 
neutralizing antibodies against RSV F”®. 

The spatial proximity of domains A and B in the S trimer allows ration- 
alization of their alternative use among coronaviruses to interact with 
host receptors. MHV uses the viral membrane distal loops decorating 
domain A to interact with CEACAM 1a (ref. 13), whereas MERS-CoV 
and SARS-CoV rely on the 3-motif protruding from domain B to bind 
to DPP4 (ref. 11) or ACE2 (refs 12 and 14), respectively (Extended 
Data Fig. 5a—d). The poor sequence conservation of the B domain 
3-motif among coronavirus S proteins, its considerable length variation 
among MHV strains (Extended Data Fig. 9) and our density-guided 
homology model of this motif indicate structural and functional differ- 
ences. These structural variations constitute the molecular basis under- 
lying coronavirus species specificity and cell tropism using a single S 
architectural scaffold. 

Sequence comparisons indicate that the MHV spike S; and S, sub- 
units respectively share ~25% and ~40% sequence similarity with 
many other coronavirus S proteins (Extended Data Fig. 9). Therefore, 
the structure reported here is representative of the architecture of other 
coronavirus S such as those of MERS-CoV and SARS-CoV. This hypoth- 
esis is further supported by the structural similarity of (1) the MHV” 
and bovine coronavirus!’ A domains; (2) the MHV, MERS-CoV"!," 
SARS-CoV" and HKU4 (ref. 29) B domains (Extended Data Fig. 10); 
(3) the post-fusion cores of MHV’, SARS-CoV®"" and MERS-CoV’; and 
(4) the isolation of infectious coronaviruses featuring a deletion of the 
A domain and using domain B as the receptor-binding domain*’. Our 
results now provide a framework to understand coronavirus entry and 
suggest ways for preventing or treating future coronavirus outbreaks. 


Online Content Methods, along with any additional Extended Data display items and 
Source Data, are available in the online version of the paper; references unique to 
these sections appear only in the online paper. 
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METHODS 


No statistical methods were used to predetermine sample size. 
Plasmids. A human codon-optimized gene encoding the MHV spike gene 
(UniProt: P11224) was synthesized with an Arg717Ser amino acid mutation 
to abolish the furin cleavage site at the S;—S, junction (S, cleavage site). From 
this gene, the fragment encoding the MHV ectodomain (residues 15-1231) was 
PCR-amplified and ligated to a gene fragment encoding a GCN4 trimerization 
motif (IKRMKQIEDKIEEIESKQKKIENEIARIKKIK)?“4, a thrombin cleavage 
site (LVPRGSLE), an 8-residue long Strep-Tag (WSHPQFEK) and a stop codon. 
This construct results in fusing the GCN4 trimerization motif in register with the 
HR2 helix at the C-terminal end of the MHV S-encoding sequence. This gene was 
cloned into the pMT/BiP/V5/His expression vector (Invitrogen) in frame with 
the Drosophila BiP secretion signal downstream the metallothionein promoter. 
The D1 domain of mouse CEACAM 1a (residues 35-142; gb NP_001034274.1) 
was amplified by PCR and cloned into a mammalian expression plasmid”, in 
frame with a CD5 signal sequence at the 5’ end, and with a sequence encoding a 
thrombin cleavage site, a glycine linker and the Fc domain of human IgG] at the 
3’ end, creating the pCD5-MHVR-T-Fc vector. 
Production of recombinant CEACAM 1a ectodomain by transient transfection. 
293-F cells were grown in suspension using FreeStyle 293 Expression Medium 
(Life technologies) at 37°C in a humidified 5% CO, incubator on a Celltron shaker 
platform (Infors HT) rotating at 130r.p.m. (for 1 | culture flasks). Twenty-four 
hours before transfection, cell density was adjusted at 1.5 x 10° cells ml~’, and 
culture grown overnight in the same conditions as mentioned above to reach 
~2.5 x 10° cells ml! the day of transfection. Cells were collected by centrifugation 
at 1,250 r.p.m. for 5 min, and resuspended in fresh FreeStyle 293 Expression 
Medium (Life technologies) without antibiotics at a density of 2.5 x 10° cells ml“. 

To produce recombinant CEACAM 1a ectodomain, 400 .g of pCD5-MHVR- 
T-Fc vector (purified using EndoFree plasmid kit from Qiagen) were added to 
200 ml of suspension cells. The cultures were swirled for 5 min on shaker in the 
culture incubator before adding 91g ml! of Linear polyethylenimine (PEI) solu- 
tion (25 kDa, Polysciences). Twenty-four hours after transfection, cells were diluted 
1:1 with FreeStyle 293 Expression Medium and the transfected cells were cultivated 
for 6 days. Clarified cell supernatants were concentrated tenfold using Vivaflow 
tangential filtration cassettes (Sartorius, 10-kDa cut-off) before affinity purification 
using a Protein A column (GE LifeSciences) followed by gel filtration chromatog- 
raphy using a Superdex 200 10/300 GL column (GE Life Sciences) equilibrated in 
20mM Tris-HCl, pH 7.5, 100 mM NaCl. The Fc tag was removed by trypsin cleav- 
age in a reaction mixture containing 7 mg of recombinant CEACAM 1a ectodomain 
and 5g of trypsin in 100 mM Tris-HCl, pH 8.0 and 20mM CaCl). The reaction 
mixture was incubated at 25°C overnight and re-loaded in a Protein A column to 
remove uncleaved protein and the Fc tag. The cleaved protein was further purified 
by gel filtration using a Superdex 75 column 10/300 GL (GE Life Sciences) equili- 
brated in 20 mM Tris-HCl, pH 7.5, 100 mM NaCl. The purified protein was quan- 
tified using absorption at 280 nm and concentrated to approximately 10 mg ml". 
Production of recombinant MHV S ectodomain in Drosophila S2 cells. To 
generate a stable Drosophila S2 cell line expressing recombinant MHV spike ecto- 
domain, we used Effectene (Qiagen) and 2 jg of the plasmid encoding the MHV 
S protein ectodomain. A second plasmid, encoding blasticidin S deaminase was 
cotransfected as dominant selectable marker. Stable MHV S ectodomain expressing 
cell lines were selected by addition of 101g ml“! blasticidin S (Invivogen) to the 
culture medium 48 h after transfection. 

For large-scale production of MHV S ectodomain the cells were cultured 
in spinner flasks and induced by 5 4M CdCl, at a density of approximately 10’ 
cells per ml. After a week at 28 °C, clarified cell supernatants were concentrated 
40-fold using Vivaflow tangential filtration cassettes (Sartorius, 10-kDa cut-off) 
and adjusted to pH 8.0, before affinity purification using StrepTactin Superflow 
column (IBA) followed by gel filtration chromatography using Superose 6 10/300 
GL column (GE Life Sciences) equilibrated in 20mM Tris-HCl, pH 7.5, 100 mM 
NaCl. The purified protein was quantified using absorption at 280 nm and con- 
centrated to approximately 4mg ml7’. 
SEC-MAILS. For size exclusion chromatography coupled with multi-angle light 
scattering (SEC-MALS) analysis, samples (0.2 ml at 1 mg ml’) were loaded onto 
a Superdex 200 10/300 GL column (GE Life Sciences, 0.4m] min“ in gel filtration 
buffer) and passed through a Wyatt DAWN Heleos II EOS 18-angle laser photom- 
eter coupled to a Wyatt Optilab TrEX differential refractive index detector. Data 
were analysed using Astra 6 software (Wyatt Technology Corp). 
MicroScale Thermophoresis. Solution MicroScale Thermophoresis (MST) 
binding studies were performed using standard protocols on a Monolith NT.115 
(Nanotemper Technologies). In brief, recombinant CEACAM 1a ectodomain 
protein was labelled using the RED-NHS (Amine Reactive) Protein Labelling 
Kit (Nanotemper Technologies). The MHV S ectodomain protein was serially 
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diluted in 20mM Tris-HCl, pH 7.5, 100 mM NaCl and the labelled recombinant 
CEACAM 1a was added to a final concentration of 500 nM before overnight 
incubation at 4°C. The CEACAM 1a concentration was chosen such that the 
observed fluorescence was approximately 1,000 U at 40% LED power. The samples 
were loaded into standard-treated Monolith capillaries and were measured by 
standard protocols using a Monolith NT.115, NanoTemper. The changes in the 
fluorescent thermophoresis signal were plotted against the concentration of 
the serially diluted MHV spike protein, and Kg values were determined using the 
NanoTemper analysis software. 

CryoEM sample preparation and data collection. Three microlitres of MHV 
spike at 1.85 mg ml“! was applied to a 1.2/1.3 C-flat grid (Protochips), which had 
been glow-discharged for 30s at 20 mA. Thereafter, grids were plunge-frozen in 
liquid ethane using a Gatan CP3 and a blotting time of 3.5s. Data were acquired 
using an FEI Titan Krios transmission electron microscope operated at 300 kV 
and equipped with a Gatan K2 Summit direct detector. Coma-free alignment was 
performed using the Leginon software*’. Automated data collection was carried out 
using Leginon™ to control both the FEI Titan Krios (used in microprobe mode at a 
nominal magnification of 22,500 x ) and the Gatan K2 Summit operated in counted 
mode (pixel size: 1.315 A) at a dose rate of ~9 counts per physical pixel per s, which 
corresponds to ~12 electrons per physical pixels per s (when accounting for coinci- 
dence loss*°). Each video had a total accumulated exposure of 53e A * fractionated 
in 38 frames of 200 ms (yielding movies of 7.6 s). A data set of ~1,600 micrographs 
was acquired in a single session using a defocus range of between 2.0 and 5.0,1m. 
CryoEM data processing. Whole-frame alignment was carried out using the soft- 
ware developed previously*°, which is integrated into the Appion pipeline”®, to 
account for stage drift and beam-induced motion. The parameters of the micro- 
scope contrast transfer function were estimated for each micrograph using ctffind3 
(ref. 37). Micrographs were manually masked using Appion to exclude the visible 
carbon supporting film for further processing. Particles were automatically picked 
in a reference-free manner using DogPicker*®. Extraction of particle images was 
performed using Relion 1.4 with a box size of 320 pixels” and applying a windowing 
operation in Fourier space to yield a final box size of 288 pixels” (corresponding 
to a pixel size of 1.46 A). From the 1.2 million particles initially picked, a subset of 
50,000 particles were randomly selected to generate class averages using RELION”’. 
An initial 3D model was generated using OPTIMOD” within the Appion pipeline. 
The entire data set was subjected to 2D alignment and clustering using RELION 
and particles belonging to the best-defined class averages were retained (~500,000 
particles). These ~500,000 particles were then subjected to RELION 3D classi- 
fication with four classes (using cl symmetry) starting with our initial model 
low-pass filtered to 40 A resolution. We subsequently used the ~230,000 best 
particles (selected from the 3D classification) and the map corresponding to the 
best 3D class (low-pass filtered at 40 A resolution) to run Relion 3D auto-refine (c3 
symmetry), which led to a reconstruction at 4.4 A resolution. We used the particle 
polishing procedure in RELION 1.4 to correct for individual particle movement 
and radiation damage*!*”. A second round of 3D classification with 6 classes 
(c3 symmetry) was performed using the polished particles resulting in the selection 
of 82,000 particles. A new 3D auto-refine run (c3 symmetry) using the selected 
82,000 particles and the map corresponding to the best 3D class (low-pass filtered 
at 40 A resolution) yielded a map at 4.0 A resolution following post-processing in 
RELION. The final map was sharpened with an empirically determined B factor 
of —220 A? using Relion post processing. Reported resolutions are based on the 
gold-standard Fourier shell correlation (FSC) = 0.143 criterion®’, and Fourier shell 
correction curves were corrected for the effects of soft masking by high-resolution 
noise substitution®*. The soft mask used for FSC calculation had a 10 pixel cosine 
edge fall-off. The overall shape and dimensions of our reconstruction agree 
with previous data although the HR2 stem connecting to the membrane is not 
resolved’”. 

Model building and analysis. Fitting of atomic models into cryoEM maps was 
performed using UCSF Chimera*® and Coot'*“°. We initially docked the MHV 
domain A structure (PDB 3R4D) and used a crystal structure of a bovine corona- 
virus domain A (PDB 4H14) to model the three-stranded 3-sheet and the a-helix 
present on the viral membrane proximal side of the galectin-like domain. Next, the 
MERS-CoV domain B crystal structure (PDB 4KQZ) was also fit into the density, 
and rebuilt and refined using RosettaCM*’. Although we could accurately align 
the sequences corresponding to the core B-sheet of the MHV and MERS-CoV B 
domains, the ~100 residues forming the 8-motif extension (residues 453-535, 
MERS-CoV/SARS-CoV receptor-binding moiety) could not be aligned with con- 
fidence. We used RosettaCM to build models of each of the 945 possible disulfide 
patterns into the density for domain B. For each disulfide arrangement, 50 models 
were generated, and there was a very clear energy signal for a single such arrange- 
ment (Extended Data Fig. 3k). Then, 1,000 models with this disulfide arrange- 
ment were sampled, and the lowest energy model (using the Rosetta force field 
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augmented with a fit-to-density score term) was selected. Owing to the poor quality 
of the reconstruction at the apex of the S trimer, the confidence of the model is 
lowest for the segment corresponding to residues 453-535, as homology modelling 
was used to fill in details missing in the map. 

A backbone model was then manually built for the rest of the S polypeptide 
using Coot. Sequence register was assigned by visual inspection where side chain 
density was clearly visible. This initial hand built model was used as an initial 
model for Rosetta de novo’. The Rosetta-derived model largely agreed with the 
hand-built model. Rosetta de novo successfully identified fragments allowing to 
anchor the sequence register for domains C and D as well as for helices a21-a95. 
Given these anchoring positions, RosettaCM*’ augmented with a novel density- 
guided model-growing protocol was able to rebuild domains C and D in full. The 
final model was refined by applying strict non-crystallographic symmetry con- 
straints using Rosetta’®. Model refinement was performed using a training map 
corresponding to one of the two maps generated by the gold-standard refinement 
procedure in Relion. The second map (testing map) was used only for calculation 
of the FSC compared to the atomic model and preventing overfitting*®. The quality 
of the final model was analysed with Molprobity””. Structure analysis was assisted 
by the PISA®® and DALI”! servers. The sequence alignment was generated using 
MultAlin°? and coloured with ESPript*’. All figures were generated with UCSF 
Chimera*. 
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Extended Data Figure 1 | Biophysical characterization of the MHV molecular mass (expressed in Da, left ordinate axis). b, MHV S binds with 
S ectodomain. a, The MHV S molecular mass was determined to be high-affinity to the soluble mouse CEACAM 1a receptor. Thermophoresis 
463.2 + 0.3 kDa (mean + s.e.m.) (corresponding to a trimer) using signal plotted against the MHV S concentration. The dissociation constant 
size-exclusion chromatography coupled in-line with multi-angle light (Ka) was determined to be 48.5 + 3.8nM. Values correspond to the average 
scattering and refractometry. The blue line represents the normalized of two independent experiments. The concentration of CLACAM 1a used 
refractive index (right ordinate axis) and the red line shows the estimated was 500 nM. 
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Extended Data Figure 2 | CryoEM analysis of the MHV S trimer. model/map (red) Fourier shell correlation (FSC) curves. The resolution 
a, b, Representative electron micrograph (defocus: 4.6 1m) (a) and class was determined to 4.0 A. The 0.143 and 0.5 cut-off values are indicated by 
averages (b) of the MHV S trimer embedded in vitreous ice. Scale bars: horizontal grey bars. 


573A (micrograph) and 44 A (class averages). c, Gold-standard (blue) and 
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Extended Data Figure 3 | CryoEM density for selected regions of the protomer (j), coloured according to local resolution determined with the 
MHV S reconstruction, local resolution analysis and density-guided software Resmap. We interpret Resmap results as a qualitative (rather than 
homology modelling of residues 453-535. The atomic model is shown quantitative) estimate of map quality. k, Rebuilding of the MHV S domain 
with the corresponding region of the map. a, b, Upstream helix. B using RosettaCM. Plot showing the energy mean and s.d. of the models 
c-e, Helix belonging to domain A (residues 284-296). f-h, Core 3-sheet. corresponding to the 30 lowest energy disulfide arrangements (out of 945) 


i, j, CryoEM density corresponding to the MHV S trimer (i) and a single for domain B. 
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Data collection 


Number of particles 82,000 

Pixel size (A 1.315 (rescaled to 1.46 
Defocus range (um 2-5 

Voltage (kV) 300 
Electron dose (e/A 53 

Refinement 
Resolution 4.0 
Map sharpening B factor (A?) -220 


Model validation 


Molprobity score (percentile) 150 (94) 
All-atom clashscore (percentile) 3.68 (97) 


Poor rotamers (%) 0 
Ramachandran favored (%) 94.26 
Ramachandran allowed (%) 99.44 
Ramachandran outliers (%) 0.56 

r.m.s.d bonds (A) 0.018 

r.m.s.d angles (°) 1.8 


Extended Data Figure 4 | Refinement and model statistics. 
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Extended Data Figure 5 | Structural organization of the S, subunit. weak to allow tracing of this segment (residues 453-535), which has been 
a, Ribbon diagram showing a single S; protomer. b, Close-up view of the traced by density-guided homology modelling using Rosetta. c, d, Ribbon 
MHV S domain B. The structural motif used as a receptor-interacting diagrams of the S; trimer viewed from the side (c) and from the top 
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Extended Data Figure 6 | Mechanisms of membrane fusion promoted Q23 as a single extended helix. c, Ribbon diagram of the SARS-CoV post- 
by coronavirus S glycoproteins. a, Ribbon diagram of the MHV S; pre- fusion HR1 helix obtained by X-ray crystallography (PDB 1WYY). The 
fusion structure. Disulfide bonds are shown as green sticks. b, Topology residue numbers corresponding to the MHV A59 sequence are indicated. 
diagram of the MHV S, pre-fusion structure. PP, di-proline that willactas —_d, Topology diagram showing the expected coronavirus S post-fusion 
a helix breaker. The presence of these di-proline motifs indicates that the conformation derived from our MHV S structure and the SARS-CoV 
post-fusion HR1 coiled-coil could not extend up to the fusion peptide as post-fusion core crystal structure shown in c. e, Ribbon diagram of a 
a single helix. This hypothesis is further supported by the observation of model of the MHV S, post-fusion conformation. Residues belonging 
a conserved disulfide bond formed between residues Cys894 and Cys905 tO Q21, O22, O23, B4g, O24 and Qy5 are not represented owing to a lack of 
(labelled 14 in a and b), which will prevent refolding of helices a2 and structural information. 
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Extended Data Figure 7 | Structural organization of the S, fusion machinery. a, Ribbon diagram of the trimer of central helices. b, c, Ribbon diagrams 
of the S, trimer (starting at residue 755) viewed from the side (b) and from the bottom (looking towards the host cell membrane) (c). Residues Ala994 
and Leu1062, which are discussed in the text, are shown in stick format. 


© 2016 Macmillan Publishers Limited. All rights reserved 


RESEARCH @aSaia 


HIV Env 


Y opis x rs 


>) 
, i) 
y, > vA Ge: 
(== M NY) (2_y) 
WS < y Z PE 


/ y 
ee 


SWES IY SS 
a xT 2 —_ 
a> 


cnt? 


fg 
Wil 
4 a > 


WW 
(W—= fae) AS 
Nees X <p} ue Se 
‘\ ~aey PAY AI IN 
WSS Sere 
\ IS SS c - fi 


sy, 
<S) 


f Y/ 


"e 
Xa 


4 
ay 


* 


Extended Data Figure 8 | Class I viral fusion proteins with exposed fusion peptide. a, MHV S (residues 870-887). b, Parainfluenza virus 5 F (PIV5 F, 
residues 103-128, PDB 2B9B). c, HIV-1 gp41 (residues 518-528, PDB 4T VP). The trimeric fusion proteins are shown as grey ribbon diagrams with the 


fusion peptides rendered in magenta. 
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Extended Data Figure 9 | See next page for eae caption. 
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Extended Data Figure 9 | Sequence conservation among coronavirus S 
glycoproteins. a, Sequence alignment of coronavirus S proteins. Bovine- 
CoV, bovine respiratory coronavirus AH187 (gi 253756585); HKU1, 
human coronavirus HKU1 (gi 545299280); HKU4, tylonycteris bat 
coronavirus HKU4 (gi 126030114); HKUS5, pipistrellus bat coronavirus 
HKUS (gi 126030124); MERS-CoV, Middle East respiratory syndrome 
coronavirus (gi 836600681); MHV-A59, mouse hepatitis virus A59 

(gi 1352862); MHV-JHM, mouse hepatitis virus JHM (gi 60115395); 
MHV-2, mouse hepatitis virus 2 (gi 5565844); OC43, human coronavirus 
OC43 (gi 744516696); SARS-CoV, severe acute respiratory syndrome 
coronavirus ZJ01 (gi 39980889); Waterbuck-CoV, waterbuck coronavirus 


US/OH-WD358-TC/1994 (gi 215478096). Asparagine residues featuring 
N-linked glycan chains visible in the MHV S reconstruction are indicated 
with a star. The S2 and S,’ cleavage sites are indicated with scissors at 
positions corresponding to the MHV S sequence. Cysteine residues 
involved in the formation of disulfide bonds are numbered according to 
Supplementary Table 2. The secondary structure elements observed in our 
MHYV S reconstruction are indicated above the sequence. The black dotted 
lines above the sequence indicate regions poorly defined in the density. 
Although the viral membrane distal loops of the A domains are weakly 
defined in the density, the availability of a crystal structure of this domain 
from the same virus (PDB 3R4D) helped with the modelling. 
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Extended Data Figure 10 | Structural similarity of B domains among coronavirus S glycoproteins. a, MHV (pink). b, MERS-CoV (orange, PDB 
4KQZ). c, SARS-CoV (red, PDB 2AJF). d, HKU4 (blue, PDB 4QZV). 
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